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ABSTRACT 

Unlike at lower redshift, where there is a 40% detection rate, surveys for 21 -cm absorption 
aris ing within the hos ts of z > 1 radio galaxies and quasars have been remarkably unsuccess- 
ful. ICurran et al.l (12008 ) suggest that this is due to the high redshift selection biasing towards 
the most optically bright objects (those most luminous in the ultra-violet in the rest-frame), 
where the gas is ionised by the active galactic nucleus. They therefore argue that there must 
be a population of fainter objects in which the hydrogen is not ionised and which exhibit a 
similar detection rate as at lower redshifts. In order to find this "missing" gas at high redshift, 
we have therefore undertaken a survey of z > 2 radio sources, selected by optical faintness. 
Despite having optical magnitudes which indicate that the targets have ultra-violet luminosi- 
ties below the threshold where all of the gas is ionised, there were no detections in any of the 
eight sources for which useable data were obtained. Upon an analysis of the spectral energy 
distributions, ionising photon rates can only be determined for three of these, all of which 
suggest that the objects are above the highest luminosity of a current 21 -cm detection. The 
possibility that the other five could be located at lower photon rates cannot be ruled out, al- 
though zero detections out of five is not statistically significant. Another possible cause of the 
non-detections is that our selection biases the sample towards sources which are very steep 
in the radio band, with a mean spectral index of (a) = —1.0, cf. —0.3 for both the 21 -cm 
detections and UV luminous non-detections. This adds the further possibility that the sources 
have very extended emission, which would have the effect of reducing the coverage by the 
putative absorbing gas, thus decreasing the sensitivity of the observation. 

Key words: galaxies: active - quasars: absorption lines - radio lines: galaxies - ultra violet: 
galaxies - galaxies: stellar content - galaxies: high redshift 



1 INTRODUCTION AND SAMPLE SELECTION 

Redshifted Hi 21 -cm provides a probe of the most abundant ele- 
ment in the Universe, through surveys which are not subject to the 
flux limitations of optical surveys. In absorption, 21 -cm probes the 
cool component of the neutral gas, which is the reservoir of raw ma- 
terial ultimately responsible for all star formation. Furthermore, the 
strength of the absorption depends only upon the column density 
of the absorber and the flux of the background source, making the 
21-cm transition potentially detectable up to redshifts of 2 ~ 50, 
where the ionosphere begins to affect the < 30 MHz radio waves. 

Despite this, redshifted 21-cm absorption is currently rare, 
with only 78 absorbers known at z ^ 0.1, 80% of which are de- 
tected at z 1. Much of this is due to the past availability of inter- 
ference free bands at low frequencies, although there are additional 
selection effects at play: 
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- For absorbers intervening the sight-lines to more distant ra- 
dio sources, the 60% detection rate at 2: ^ 1, compared to 20% 
at 2 ^ 1, in systems known to have high hydrogen column den- 
sitiefl can be attributed to the geometry effects introduced by an 
expanding Universe, where the coverage of the background flux is 
systematicall y lower at redshifts of 2 ;S 1 jCurran & Webblliooel: 
ICurranll2012l) . 

- For absorbers associated with the radio source itself, the de- 

tection rate at 2 1 is 40%, compared to 17% at 2 ;S 1 

dCurran & Whitindl2010l) . This is believed to be due to the higher 
redshifts biasing towards the most ultra-violet luminous objects, 
where the intense flux ionises the gas. Specifically, associated 21- 
cm absorption has never been detected where the A = 1216 A con- 
tinuum luminosity of the active galactic nucleus (AGN) exceeds 
L1216 ~ 10^=* WHz-i l lCurran et al...2008i and recently confirmed 



^ Damped Lyman-o absoiption systems (DLAs), with neutral hydrogen 
column densities of A^hi > 2 x 10^" cm~^. 
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by iGrasha & Darlind I20T1I) . At these frequencies (2.47 x 10^^ 
Hz) and above, the photons have enough energy to excite the hy- 
dro gen beyond the ground st ate, so that it cannot absorb in 21- 
cm. ICurran & Whitind \2Q\% have extended this to the ionising 
(A ^ 912 A) radiation and show that a similar (i.e. Lqi2 ~ 10^'^ 
W Hz~^) cut-off also applies. The large fraction of non-detections 
at high redshift are therefore due to the flux limited nature of the 
optical surveys selecting the brightest objects, where the observed- 
frame optical light is rest -frame ultra-violet, in wh ich all of the gas 
is believed to be ionised jCurran & Whitindl2012h . 

In order to find the neutral gas missing from high redshift 
AGN, we suggest that high redshift 21 -cm surveys should be di- 
rected towards the most optically faint sources. We have there- 
fore embarked on an observing campaign of objects se l ected b y 
faint optical magnitudes. As discussed in ICurran et alj j2011dh . 
our usual sourc e catalogue, the Parkes Half-Jansky Flat-spectrum 
Sample (PHFS.^Drink water et alj|l997l) . yielded only two sources 
with faint blue magnitudes {B > 22) in the 90-cm band (z — 
3.09 — 3. 63il, both of whic h have been previously searched for 
in 21-cm dCurran et allbOOSl) . We therefore compiled all available 
radio catalogues, which give both redshifts and magnitudes, and 
selected those with B ,V, R ox I magnitudes which indicate that 
1-1216 ^ 10^^ W Hz"^ at the given redshift (Fig.[D. We further 
shortlisted those in which 21-cm is redshifted into the 290-395 
MHz band of the Green Bank Telescope (GBT) or the UHF-low 
band (250-460 MHz) of the Westerbork Synthesis Radio Telescope 
(WSRT) and which reached sufficiently high elevations at these lo- 
cations, while having flux densities estimated to be 0.2 Jy at the 
redshifted 21-cm frequency[f| This gave eleven sources which we 
observed (Table [TJ, eight of which were not completely ruined by 
RFI, and as seen from Fig. [T] all are believed to be below the crit- 
ical luminosity. For the « 40% detection rate at Z/1216 ~ 10^^ 
W Hz~^, we therefore expect approximately three detections, al- 
though there were none. We discuss possible reasons for the exclu- 
sive non-detections in this paper. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Green Bank observations 

Each of the sources targetted with the Green Bank Telescope were 
observed for a total of two hours on source, with the aim of 
reaching an r.m.s. noise level of J; 5 mjy per 10 kms^^ chan- 
nel (for Tsys = 70 K). This gives a 3cr optical depth limit of 
r ~ 0.03 per channel for a flux density of 0.5 Jy, i.e. a sensitiv- 
ity to iVm ~ 3 X 10^^(Tspin//) cm"^, the lower limit for most 
of the published searches (see figure 4 of ICurran et alj[20o3) . We 
used the Rcvr_3 4 2 receiver backed with the GBT Spectrometer in 
narrow bandwidth, high resolution mode in 4 x 12.5 MHz bands, 
while maintaining a redshift coverage of Az ~ ±0.06 in each 
band. Each source was observed in two orthogonal linear polari- 
sations which, after flagging of bad data, were averaged together 
(using the GBTIDL software): 

NVSS J012142+132058 was observed at 314.53 MHz on 17 June 
2010, over 4096 channels, giving a spacing of 3.052 kHz (2.92 

Applying the B — z curve (figure 5 of lCurran et alj200g|) . gives B C 22 
for L1216 i 10^3 W Hz-l at these redshifts. 

Except for J06 17+50 12 for which no fluxes were available from the 
NASA/IPAC Extragalactic Database, although we report 1.28 Jy at 342 
MHz here. 



km s~^). After flagging of the worst RFI affected scans, 0.6 hours 
of data remained with an average system temperature of Tsys = 90 
K, although some spikes were still present in the bandpass (Fig.|2]l. 
NVSS J020510+224250 was observed at 315.23 MHz on 31 Au- 
gust 2010, over 4096 channels, giving a spacing of 3.052 kHz (2.91 
km s^^). Tsys = 87 K and the YY polarisation exhibited negative 
fluxes for the whole observations and after flagging out the worst 
RFI affected on the XX polarisation 0.68 hours of data remained. 
NVSSJ023111+360027 was observed at 348.22 MHz on 1 
September 2010, over 8196 channels, giving a spacing of 1.526 
kHz (1.31 kms~^). After flagging, 0.74 hours of data remained, 
although all of the scans exhibited negative fluxes. Tsys = 75 K 
and between the RFI spikes, the r.m.s. noise levels were low (4 
mJy per 10 km s~^), although the spikes raised the mean to 189 
mJy over the observed range. 

WNJ0617+5012 was observed at 342.02 MHz on 6 July 2010, 
over 8196 channels, giving a spacing of 1.526 kHz (1.34 kms"^). 
The mean system temperature was 76 K and after flagging, 0.63 
hours of data remained. There is an "absorption feature" close to 
the expected redshift of 2; = 3.153 (fitted by two narrow Gaus- 
sians at an offset of 220 and 241 km from km s^^, defined 
by 342.02 MHz), which is present on all four sub-averages of the 
data. However, this is also true for the two "emission features" at 
~ —1500 and —1100 kms~^ and so we must attribute these to 
RFI. Again, the regions between the spikes had low r.m.s. noise 
levels (5 mJy per 10 kms~^), although in Table [T] we quote the 
value over 337.56 -343.17 MHz, were the spikes are fairly weak. 
PKS 0742+10 was observed at 390.30 MHz on 13 August 2010, 
over 16 384 channels, giving a spacing of 763 kHz (0.59 kms^^). 
After flagging of the worst RFI affected scans, 0.6 hours of data re- 
mained with an average system temperature of Tsys = 81 K. How- 
ever, intermittent spikes remained, meaning that the r.m.s. noise 
level quoted in Table [T] only applies to the clean regions between 
these, with a spike being present over z = 2.6302 — 2.6309 
(391.20-391.27 MHz), very close to the emission redshift of the 
source. 

FIRST J074729.3+365438 was observed at 355.81 MHz on 7 June 
2010, over 16 384 channels, giving a spacing of 763 kHz (0.66 
km s~^). The mean system temperature was 63 K and 0.8 hours 
of data remained after flagging the integrations affected most by 
RFI. As seen in Fig.|2] the bandpass was still dominated by spikes 
and so in Table[T]we quote the r.m.s. noise level over the relatively 
clear 355.78-357.15 MHz. 

NVSSJ092022-071217 was observed at 377.77 MHz on 9 
September and 29 December 2010, over 16 384 channels, giving a 
spacing of 763 kHz (0.61 km s~^). The mean system temperature 
was 65 K and flagging of badly affected data gave a total integration 
time of 0.63 hours. Even so, the bandpass was still dominated by 
spikes with RFI close to the expected frequency of the redshifted 
21-cm line preventing us from assigning a limit. 
B2 1121+31B was observed at 336.83 MHz in three sessions span- 
ning from 6 July to 5 August 2010. The band was split over 4096 
channels, giving a spacing of 3.052 kHz (2.73 km s"^). The mean 
system temperature was 62 K and RFI was minimal, with 0.90 
hours of data remaining after flagging the worst of this. Despite 
the absence of severe RFI and the flat bandpasses, like J0205+2242 
each scan exhibited a negative flux, in both polarisations, over the 
whole run. Negative fluxes can often caused by a region of high 
confusion leading to the off-measurement being stronger than the 
source. At 337 MHz, the half-power beam width of the GBT is 
~ 40', within which the NASA/IPAC Extragalactic Database lists 
17 000 sources, 140 of which are classified as radio sources. There- 
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Figure 1. The r ed magnit ude (top) and A = 1216 A continuum luminosity (bottom) for the redshifted associated 21-cm searches (compiled in 
ICuiTan & Whitin g 2010; All ison et al.ll20lA . The filled symbols show the detections and the unfilled symbols the non-detections, with the asterisks showing 
our targets (all of which are non-detections), i.e. the eleven points in the top panel and the eight for which useful data were obtained in the bottom panel (with 
ii2i6 as originally estimated and not corrected for extinction, as the remainin g points are). In the top panel the curve shows which B magnitude corresponds 
to L1216 = 10^'^ W Hz~^ for a spectral slope of a = —1.5 (cf. figure 5 of lCurran et aljlioo^ used to select the targets and in the bottom panel the line 
signifies L1216 = 10^'^ W Hz^^, above which associated 21-cm has never been detected. Note that for J023 1-1-3600, at 2 = 3.079, / = 25.0 is shown (Table 



fore confusion within the beam/sidelobes/off-position leading to a 
negative flux is a possibility. Although the flux levels may be not be 

S iroperly ca librated, the data themselves are believed to the reliable 
Riesll2012l) . as is apparent from the consistency between scans for 
both J0205-H2242 and 112 1-1-3 IB. 

2.2 Westerbork observations 

For the WSRT observations we requested four hours for each 
source using all 15 antennas (giving 105 baseline pairs), with 
the aim of reaching an r.m.s. noise level of ^ 20 mjy per 10 
kms~^ channel. For a flux density of 0.5 Jy, this gives a So- 
optical depth limit of r ~ 0.1 per channel, a sensitivity to 
A^Hi ^ 2 X 10^*(rspin//) cm"^ We used the UHF-low receiver 
backed with the correlator over 2x5 MHz bands, each in two 
polarisations over 1024 channels, in order to have a channel 
spacing of ~ 4 km s^^, while maintaining a redshift coverage of 
Az ~ ±0.04. For each source the delays were all self calibrated 
with 3C 48, 3C 147 and 3C 286 being used for bandpass calibration 
and the data were reduced using the MIRIAD interferometry 
reduction package, with a spectrum being extracted from each 
cube. As per the GBT spectra, these are shown in Fig. [2] and 
summarised in Table[T] 

NVSSJ012142-I-132058 was observed on 5 January 2010 for a 



total of 5.0 hours. Severe RFI meant that extensive flagging of the 
data was required, leaving 45 baseline pairs. Despite the flagging, 
a reasonable image could not be produced and so a spectrum 
was obtained by averaging together the remaining baselines pairs, 
giving a much reduced flux with two spikes remaining. An r.m.s. 
noise level of 10 mJy per 10 km s^^ channel was reached between 
these spikes, which rises to 134 mJy averaged over the whole 
bandpass. Each unsmoothed channel was 4.66 kms~^ wide and 
the synthesised beam was 772" x 83" . 

NVSSJ020510-I-224250 was observed on 19 February 2010 
for a total of 5.0 hours. Again severe RFI required extensive 
flagging, leaving only the XX polarisation and 78 baseline pairs 
which produced a poor quality image. A spectrum was, however, 
extracted from the cube, which had an r.m.s. noise level of 77 mJy 
per 10 km s^^ channel between the two main spikes, which rises 
to 130 mJy over the whole bandpass. Each unsmoothed channel 
was 4.64 km s^^ wide and the synthesised beam 237" x 61". 
NVSS J023111-I-360027 was observed on 10 February 2010 for a 
total of 3.5 hours. Again, RFI required extensive flagging, leaving 
only 55 baseline pairs. From this an image of the target could not 
be produced, with patchy emission over the field being evident as 
a ripple in the extracted spectrum. Each unsmoothed channel was 
4.20 km s^^ wide and the synthesised beam 256" x 69". 
B20300-I-37A was observed for 2.7 hours on 14 December 2009 
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for 4.5 hours. After flagging 66 baselines remained, but RFI was 
still too severe to allow us to produce an image and so we show the 
average of the unflagged baseline pairs. Each unsmoothed channel 
was 3.61 km wide. 

WN J0617+5012 was observed at 342.02 MHz for 6.0 hours on 18 
January 2010, but severe RFI meant that no data could be used. 
PKS 0742+10 was observed at 390.30 MHz for 5.8 hours on 
14 December 2009, but again, severe RFI marred the whole 
observation. 

FIRST J074729.3+365438 was observed at 355.81 MHz for 6.0 
hours on 1 1 March 2010. After flagging the worst RFI, 55 baseline 
pairs remained from which an image of the target could not be 
produced. Nevertheless, the spectrum extracted from the field was 
relatively spike free and, as in the case of J023 1+3600, there is a 
strong ripple evident, although the spectrum is very noisy. Each 
unsmoothed channel was 4.10 kms^^ wide and the synthesised 
beam was 244" x 66". 

VLSS J1115.1+5016 was observed at 401.24 MHz on 12 Decem- 
ber 2009 for 4.0 hours. However, severe RFI in this band meant 
that none of the data where useable. 

B2 1121+31B was observed at 336.83 MHz for 5.2 hours on 19 
January 2010. After flagging, 66 baseline pairs remained from 
which an image was produced. The image was of poor quality and 
the extracted spectrum is over the central ± ~ 15', with the flux 
density showing a strong dependence on the extent of the image 
used, thus limiting the reliability of this. Each unsmoothed channel 
was 4.35 km wide and the synthesised beam was 288" x 62" . 
B2 1240+39 was observed at 453.66 MHz for 3.0 hours on 5 April 
2010. After flagging, 55 baseline pairs remained and a spectrum 
was extracted from the cube. Each unsmoothed channel was 3.27 
km s^^ wide and the synthesised beam was 303" x 73". 



3 POSSIBLE REASONS FOR THE NON-DETECTION OF 
21-CM 

3.1 Sensitivity limits 

Despite our deliberate selection of optically faint sources, with 
continuum luminosities below the threshold of Z/i2i6 = 10^^ 
W Hz^^, we have not detected Hi 21 -cm in any of the eight 
sources for which there w ere useable data. Below this critical value, 
ICurran & Whilh^ | |20IO|) find a ^ 40% general detection rate, 
which consists of « 50% for compact objectjf] and ~ 40% for 
others. That is, we may expect approximately three detections out 
of the eight sources for which the data were not overwhelmed by 
RFI. 

One possible reason for the non-detections is that all of the 
sources where affected to some degree by RFI, giving some rela- 
tively poor sensitivity limits (Table [TJ, where detections of optical 
dep ths of Tobs ~ 0.1 have been documen ted to be in the minor- 
ity jVermeulen et aljr2003l . lMorgantill2004h . In order to investigate 
this, in Fig. [3] we show the obs optical depths for all the redshifted 
21-cm searches. Since many of our sources were observed by both 

^ Gigahertz peaked spectrum (GPS), compact steep spectrum sources 
(CSS) and compact symmetric objects (CSO). Although the 21-cm detec- 
tion rate is generally believed to be higher for these than the general pop- 
ulation, when the Z/1216 ^ 10^^ W Hz~^ sources are removed the rates 
ai'e similar, indicating that the higher rates may be due to the generally lo w 
ultra-violet luminosities of the compact objects jCurran & Whitindl2O10l) . 



the GBT and WSRT, we use the best limit in each case and smooth 
these, as well as rebinning the 3(t limits from the literature, to a 
channel spacing of 167 kms"\ the mean FWHM of the 21-cm 
detections. This normalises all of the limits, published at various 
spectral resolutions, and gives the limit for the detection of putative 
line of width FWHM = 167 km s^^ within a single channel. 

This normalisation has the effect of moving all of our limits to 
Tobs ~ 0.1, while degrading the column density limits (see Fig.O 
by a factor of 4 to those quoted in Table [T] From Fig. [3] we see 
that our limits are located in the densest clustering of the detections 
and so for the rest of the analysis, since we are attempting to explain 
the exclusive non-detections, we assume that (at least some of) the 
targets have been searched sufficiently deeply. 

We note also a weak correlation between the obs optical depth 
and the profile width (left panel]fl which may suggest that weak, 
wide profiles may be missed by 21-cm surveys, despite having sim- 
ilar column densities to the detections (Allison et al., in prep.). For 
completeness, we also show how the optical depth is distributed 
with redshift (middle panel) and the rest frame 1 .4 GHz flux den- 
sity, from which we see the expected anti-correlation. 

3.2 Radio properties 

A possible effect which could reduce the probability of detect- 
ing 21-cm absorption is the non-deliberate selection of these 
high redshift radio galaxies having steep radio spectra (all, apart 
from B0742+I0 bein g "Ultra-Steep Spectrum Radio Sources", 
iDe Breuck eTal]|2002h . This may have an effect on the actual op- 
tical depth of the line, r, which is related to the observed opti- 
cal depth, Tobs = AS/S, through r = -In ^1 - where 
/, the covering factor, quantifies how much of the observed flux 
is intercepted by the absorber. In the optically thin regime (where 
AS ^ 0.3 S), this expression reduces to r ~ AS/{Sf) = Tohs/ f ■ 
Since, by definition, / ^ 1 then the actual optical depth limit is 
T ^ Tobs El That is, the limit is affected by how effectively the ab- 
sorbing gas intercepts the background emission. 

Unfortunately, for this sample there is only high resolution 
radio imaging available f or one source - B 0742+ 10, which sub- 
tends 1.2 mas at 15 GHz l lStanghenini et al.ll200lh . At z = 2.630 
this corresponds to a linear extent of 10 pc (using a standard A 
cosmology with Ho — 71 km s~^ Mpc~^, ^matter ~ 0.27 and 
SIa = 0.73). Of the remainder of the sample, only four are in the 
area of sky currently covered by the Very Large Array's FIRST 
(Faint Images of the Radio Sky at Twenty Centimetres) survey, 
having deconvolved minor axes of 0.9I"(7 kpc) for J0747+3654, 
1.02 "(8 kpc) for J0920-07I2, I.58"(I2 kpc) for BI121+31B 
and 0.75"(6 kpc) for B 1240+39. Since these are typically un- 
resolved by the 5 "synthesised beam, assigning a source size of 
J: 40 kpc for each component may be more apt. For J0121+I320, 
J0205+2242, J023I+3600, B0300+37A and J0617+5012, which 
are all unresolved by the NRAO VLA Sky Survey (NVSS), the 
typical I7"minor axis of the beam gives source sizes of 5 130 kpc 
for these redshifts. 

Given the lack of high resolution radio imaging, information 

^ The observed distribution has a 1.4% probability of occuring by chance, 
which is significant at 2AGcr, assuming Gaussian statistics. 
® Since / is generally unknown (see Curran 20l3), we leave this (as well 
as the spin temperature) as a free parameter (Table[T). Assuming / has its 
maximum value of unity, as if often the case in the hterature, improves the 
perceived sensitivity of the observation, where in fact / < 1. 
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Figure 2. The spectra from the H I 21 -cm searches shown at a spectral resolution of 10 km . The ordinate gives the flux density [Jy] and the abscissa the 
barycentric frequency [MHz]. The downwards an'ow shows the expected frequency of the absorption from the optical redshift, with the horizontal bar showing 
a span of ±200 km s^^ for guidance. The scale along the top axis shows the redshift of H I 21-cm. 
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Figure 3. The observed optical depth versus the FWHM of the profile (left), redshift (middle) and rest-frame 1.4 GHz flux density (right). The symbols 
are as per Fig. [T] where in the left panel, for the sake of clarity, only our hmits are shown. These, and the other li mits in the remaini ng p anels, have been 
resampled to a spectral resolution of 167 km s~^, the mean of the detections which range from 18 to 475 km s"^ iCarilh et aljl 19921 and lVermeulen et alj 
I2OO3L respectively). The unfilled histogram shows the distribution for the detections and the hatched histogram the limits to our eight targets not overwhelmed 
byRFI. 



on the source size may be gleaned from the radio SEDs, given 
that the turno ver frequency is anti-correlated with the size of the 
radio source l lFantietalJll99(lh . In order to explore this possibil- 
ity, in Fig. |4l we show the spectral energy distributions (SEDs) 
of our eight targets, from which the derived turnover frequen- 
cies and spectral indices are summarised in Table (2] From this 



we see that B0742-I-10 is the only source with an indisputable 
turno ver, consistent with its observed compactness (10 pc at 15 
GHz. lStanghellini et al.l200lh . Given that none of the other sources 
exhibit a clear turnover above the lowest frequency which has been 
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Table 1. The observational results for the ^ 0.2 Jy optically faint northern sources which fall into the GBT an d WSRT 90-cm bands. With 
the exception of B0742+10 (which is from IWall & Peacock|[T98^ . all are ultra-steep spectrum sources ffrom lDe Breuck et al.|[20o3) . The 
second and third columns list the magnitudes and their values. The fouilh column gives the source redshift, followed by the observed frequency 
range [MHz] over which AS, the r.m.s. noise [mJy] reached per 10 km s~l channel, is applicable. S is the continuum flux density [Jy] and 
■'oba = ~ ~ SAS/Scont) gives the 'ia Umit to the observed optical depth of the line per 10 km channel. This is followed by the 
resulting neutral hydrogen column density, where A'jji = 1.823 X lO^* [Tspml f) /fobs d,v) [cm~^ K~l], Tgpin being the spin temperature 
and / the coveting factor (see Sect. 13.2) . Lastly we Hst the redshift range over which the limit applies. 



Source 


Mag 


Value 


z 


I'obs 


Ao 


S 


'^obs 


Tel. 




z-range 


0121+1320 




> 24 


3.516 


313.43-316.40 


23 


0.455 


< 0.15 


GBT 


< 


2.7 X 


1018 


3.4893-3.5318 








312.11-317.06 


10 


0.103 


< 0.29 


WSRT 


< 


5.3 X 


1018 


3.4799-3.5510 


J0205+2242 




> 24 


3.506 


309.93-317.60 


24 


0.251 


< 0.29 


GBT 


< 


5.3 X 


1018 


3.4723-3.5830 










367.57-372.16 


77 


0.574 


< 0.51 


WSRT 


< 


9.3 X 


1018 


2.8167-2.8643 


J023 1+3600 


/ 


25.0 


3.079 


342.68-353.46 


189 






GBT 


















345.80-350.75 


5.3 


0.377 


< 0.042 


WSRT 


< 


7.7 X 


10" 


3.0496-3.1076 


B0300+37A 


''s 


23.2 


2.506 


400.94^08.65 








WSRT 










J06 17+50 12 


R 


> 24 


3.153 


337.56-343.17 


17 


1.277 


< 0.040 


GBT 


< 


7.3 X 


lOi'^ 


3.1391-3.2079 










339.59- 344.53 








WSRT 










B0742+10 


R 


23.7 


2.630 


391.32-394.74 


22 


0.874 


< 0.076 


GBT 


< 


1.4 X 


1018 


2.5983-2.6298 










388.78-393.61 








WSRT 










J0747+3654 


■"S 


> 23 


2.992 


355.78-357.15 


15 


1.232 


< 0.037 


GBT 


< 


6.7 X 


lOi^ 


2.9771-2.9923 










353.39-358.34 


100 


0.466 


< 1.0 


WSRT 


< 


1.8 X 


1019 


2.9638-3.0194 


J0920-0712 


R 


22.4 


2.760 






0.621 




GBT 










Jl 115+5016 




> 24 


2.540 


396.21^06.21 








WSRT 










B1121+31B 




21.3 


3.217 


332.74-337.16 


7.9 




< 0.042* 


GBT 


< 


7.7 X 


lOi^ 


3.2129-3.2688 










334.56- 338.70 


66 


0.876 


< 0.23 


WSRT 


< 


4.2 X 


1018 


3.1937-3.2456 


B 1240+39 




23.6 


2.131 


451.17^56.15 


39 


0.348 


< 0.41 


WSRT 


< 


7.5 X 


1018 


2.1139-2.1483 


Notes: *Derived by 


interpolating 


! 5i5i 


MH, = 1.19 Jv JWaldram 


etalJllS 


9^ and 5365 


MHz = 


0.525 


Jv jDoualas 


etalJll996l). aivina a 



continuum flux of 5373 mHz ~ 0.57 Jy. 
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Figure 4. The SEDs of our eight targets not overwhelmed by RFI. The broken curve shows the fit to the radio data, the dotted line to the full data and the 
broken/dotted to the UV data, where possible (see the Appendix). The vertical dotted line signifies a rest-frame frequency of 3.29 X lOi^ Hz (A = 912 A). 



observed in the literatur^ll it is possible that the majority of our 
targets may have very extended radio morphologies, thus making 
T S> Tobs in the case of the limits. 

In order to verify this through the comparison of the radio 



Whether there are turnovers in the SEDs of J023 1+360, J0747+3654, 
B1121+31B and B1240+39 are disputable (Fig.|4). 



SEDs with the rest of the 21 -cm searched sources (compiled in 
ICurran & Whitinal201(]| ; [Amson et al.l2012h . we fit the radio SEDs 
(see the Appendix ) and use the ASURV survival analysis package 
jlsobe et al] 1986h to obtain mean values of the turnover frequency 
and spectral index for various sub-samples (Table (3} • From this it 
is seen that the mean turnover frequency for this sample is indeed 
lower than for any of the other sub-samples, indicating that the ra- 
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Table 2. The radio properties and luminosities of the targets, i^^^q 
is the rest-frame turnover frequency [GHz], where the upper limits 
designate no observed turnover (which is thus assumed to occur 
below the observed frequencies), followed by the spectral index at 
the rest-frame 21 -cm frequency, a. In the last three columns we 
list the A = 1216 A (as oiiginally estimated and not coiTected 
for extinction) and A = 912 A (as derived from the fits in the 
Appendix) rest-frame continuum luminosities [W Hz~^], as well 
as the ionising (A ^ 912 A) photon rate [s^^] in logs. 



Source 




a 


^1216 


^912 


rate 


J0121+1320 


< 0.331 


-1.22 


22.50 






J0205-H2242 


< 0.331 


-1.18 


22.15 






J023 1+360 


< 0.302 


-1.18 


21.28 






J06 17+50 12 


< 1.41 


-0.89 


22.00 






B0742+10 


6.46 


0.29 


22.00 


22.20 


55.68 


J0747+3654 


0.17 


-1.23 


22.54 


23.30 


57.04 


B1121+31B 


0.040 


-1.21 


22.98 


22.97 


56.03 


B 1240+39 


0.012 


-1.50 


22.60 







Table 3. The mean values of the turnover frequency [Hz] and the spectral 
index at the rest-frame 14 20 MHz for all of the sources searched in red - 
shifted 21-cm absorption jCurran & WhitindbOlOtlAlUson et alj|2012l) . 
Up gives the number of data points and n\ the number of these which are 
limits. Note that the various sub-samples from the second row on exclude 
the current sources, which are summarised in the first row. 



Sub-sample TURNOVER FREQUENCY SPECTRAL INDEX 







Tip 


ni 


(a) 


rip 


This paper 


7.91 ±0.35 


8 


4 


-1.02 ±0.18 


8 


Detections 


8.19 ±0.14 


56 


23 


-0.29 ±0.08 


58 


Non-detections 


8.03 ±0.09 


171 


92 


-0.39 ±0.04 


173 


Dets + nons 


8.07 ±0.08 


227 


115 


-0.37 ±0.04 


231 


UV luminous 


8.14 ±0.29 


19 


11 


-0.26 ±0.12 


19 



dio sources are generally larger. However, this is based upon four 
limits to the turnover frequencies out of eight SEDs and, as warned 
by ASURV, may be unreliable|f| In the case of the spectral indices, 
where there are no limits, we see that the spectra of our sample are 
significantly steeper than for the other sub-sampleslf| Thus, if the 
spectral index is a reliable tracer of the background source size, the 
fact that this sample consists almost exclusively of steep spectrum 
sources, indicating large radio source sizes, means that poor cover- 
age of the background source cannot be ruled out as the cause of the 
non-detections. Note, however, for the "UV luminous " sources (the 
19 wh ich lie above the Luv = 10^^ W Hz"^ cut-off. Icurran et al.l 
the mean spectral index is close to that of the 21-cm detec- 
tions, thus ruling out the same predisposition towards large radio 
sources being the cause of the non-detections. 



* As seen from Fig.|4] it is conceivable that seven of the turnover frequen- 
cies are limits. Furthermore, the high redshift sources will also be subject 
to a higher mean turnover due to that fact that the lowest typical value, 74 
MHz, is in the observed-frame, giving the high redshift sources higher up- 
per limits in the rest-frame. 

^ Note also that the non-detections have steeper indices than the detections, 
although these are consistent within the uncertainties. 



3.3 Ionising luminosities and photon rates 

Although all of our sources were selected to have A = 1216 A 
continuum luminosities below the critical value of I/1216 ~ 10^'^ 
W Hz~^ (Fig. [B, it is clear that the total ionising (A ^ 912 A) 
luminosity, J°° {Lv/v) dv, provides a much better measure of the 
ionising flux than a mo nochromatic value (be it I/1216 or L912, 
ICurran & Whiting|l2012h . In equilibrium, the total io nising lumi- 
nosity is related to the recombination of the atoms via (jOsterbroc^ 
Il989l) . 

poo /"^ion 

/ —^dv^4n UpUeaAr'^ dr, (1) 

^-.o„ Jo 

where h is the Planck constant, giving the number of ionising pho- 
tons per second. On the right hand side of the expression, rip and 
rie are the proton and electron densities, respectively, ua the ra- 
diative recombination rate coefficient of hydrogen and rion is the 
extent of the ionisation (the "Stromgren sphere"). For a gas which 
decreas es in density with increa sing distance from the ionising 
source, ICurran & Wliitin j ( l2012h show that non 00 for a fi- 
nite ultra-violet luminosity and that for a large spiral this value is 
L912 ^ L1216 ^ 10^'' W Hz~^, thus explaining this critical lumi- 
nosity above which 21-cm absorption has never been detected. 

From the mean SEDs of all of the associated 21-cm searches, 
the critical ionising photon rate arising from thi s luminosity is esti- 
mated to be (L^/hu) du = 2.9x10^*^ s"^ jCurran & Whitind 
I2OI2I) . From the polynomial fits (see the Appendix), we could 
determine photon rates for three of the targets presented here 
(Fig. |4j, one of which is above the estimate of the critical value 
(logio[2.9 X lO'^^] = 56.5, cf. Tabled. Applying fits to each of 
the SEDs, the highest photon rate which could be determined for 
a 21-cm detection is 5.1 x 10^^ s~^, which is below or close to 
the photon rate of the three targets (logjQ[5.1 x 10^^] — 55.7, see 
Fig. |5j- However, given that we have no estimates of the photon 
rates for the remaining five sources, from the previous estimates of 
the A — 1216 A continuum luminosities (Table[2ll, it is possible that 
several of these have rates lower than for the three for which we do 
have esti mates. For instance, J01 21 + 1320 has been detected in CO 
emission JPe Breuck et al. which may suggest a UV lumi- 

nosity below the critical value, although further studies are required 
to find if this applies to the warm emitting molecular gas. While the 
monochromatic luminosities may not provide a reliable estimate of 
the photon rates, if, for the sake of argument, the remaining five 
sources do have ionising photon rates below the critical value, for 
a 40% detection rate (Sect. 13. It , the binomial probability of zero 
out of five detections is 0.078. This is significant at 1.76ct, assum- 
ing Gaussian statistics, and so not statistically important. Recalling 
that a possible reason for the non-detections is that the targets have 
not been searched sufficiently deeply (Sect. I3.lt . this significance 
could be lower, thus not requiring us to "explain away" all eight 
non-detections. 

Lastly, given that {L^/hv)dv « 3 x lO^*^ s"^ is suffi- 
cient to ionise all of the gas in a large spiral, photon rates lower than 
this will com pletely ionise gas disks of correspondingly smaller 
scale-lengths jCurran & Whitin j2012l) . The sample presented here 
is at more than double the look-back time of the majority of 21-cm 
detections (Fig.[T) and so if there is any evolution in galactic mor- 
phology, with a larger fraction of smaller galax ies at higher redshift 
( iBaker et al.ll2000l : iLanfranchi & Friacall2003[) . we may expect the 
value of the critical photon rate to be lower at the redshifts searched 
here, making the detection of 21-cm absorption all the more diffi- 
cult. 
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Figure 5. The A = 1216 A continuum luminosity (top), A = 912 A continuum luminosity (middle) and number of ionising photons per second (bottom) 
versus the 21-cm line strength (1.823 X 10^* / rdv), see Sect. 13.31 The horizontal lines in the top and middle panels shows the critical monochromatic 
luminosity of lO'^'^ W Hz~^ and in the bottom panel JJ" {L^/hv) du = 5.1 X 10^^ s~^, the highest photon rate at which a 21-c m detection occurs 
(see main text). The symbols are as per Fig.[T] with the left pointing arrows showing the Scr upper limits to the 21-cm line strength (from lCunan & Whitind 
bOldlCurran et alj201 idjAllison et al.l20ia\ As per Fig. [3] these have been rebinne d and integrat ed over a FWHM = 167 km s-^, which has the effect of 
degrading our limits by a factor of fsi 4 in comparison to those quoted in Tableni(see lCurrai]|2012l for details). The unfilled histogram shows the distribution 
of the detection fine strengths and the hatched histogram the limits to eight targets not overwhelmed by RFI. 



Evolution may also play a role in that if the molecular (H2) 
comp onent constitutes a larger frac tion of the gas at high red- 
shift l lObreschkow & Rawlingsll2009h . lower relative column den- 
sities will make the H I more difficult to detect. However, as seen 
from Fig. [T] 21-cm has already been detected at z J: 2 and, al- 
though molecular gas has been detect ed in emissio n in over 100 
ultra-l uminous infrared galaxies (see ICurranI l2009l : ICurran et al.l 
l2011al and references therein), searches for molecular absorption 
within the hosts of m illimetre-loud sources have proved unfruitful 
( ICurranetall2011bl) PI 



Given that J0121-I-1320 is detected in CO emission and that a further five 
of the tai'gets are known exhibit strong FIR emission (Fig.|4), these sources 
should be searched in the millimetre band to rule out the possibility of a 



4 SUMMARY 

We have undertaken a survey for the hydrogen apparently miss- 
ing in high redshift radio galaxies and quasars by selecting tar- 
gets in which the optical magnitudes indicate that these are below 
the UV luminosity cut-o ff, above which all of the gas is ionised 
jCurran & Whitind2012l) . Despite this, there were no detections in 
the eight z ^ 2 objects for which useful data were obtained. Upon 

high molecular fraction suppressing the atomic abundance (if not already 
done so, although there are no published searches for molecular emission in 
the remaining targets). However, the cool absorbing molecular gas, rather 
than warm emitting co mponent, may be expected to be the coincident w ith 
the cool atomic gas (see lCombes & Wiklin3ll998l : lMurphv et alj|200ll and 
references therein). 
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an examination of the spectral energy distributions of the targets, 
we suggest two possible reasons for this: 

(i) Due to our requirement of optically faint, radio-loud objects, 
for which their redshift placed them within the GBT and WSRT 90- 
cm bands, our sample was domin ated by the ultra-steep spectrum 
sources of iDe Breuck et all l l2002h . By fitting radio SEDs to these, 
we find a mean spectral index of (a) = —1.0, which is significantly 
higher than (a) — —0.3 and —0.4 for the 21-cm detections and the 
remaining non-detections, respectively. This suggests that most of 
our targets may have extended radio emission, thus reducing the 
coverage of the background flux by the absorbing gas, lowering the 
effective sensitivity of our survey. 

(ii) Although our magnitude selection predicted that all of the 
targets are below the critical luminosity of Luv ~ 10^'^ W Hz~^, 
those for which the ionising photon rate could be determined are 
above or close to 5.1 x 10^^ ionising photons s^^, the highest value 
which can be determined for a 21-cm detection. However, the pho- 
ton rate could only be estimated for three of the eight sources and 
their monochromatic A = 1216 A continuum luminosities may 
nevertheless suggest that these are below the threshold. Without 
sufficient blue/ultra-violet photometry, this is, however, specula- 
tion. 

Although it is impossible at this time to determine which reason 
is main the culprit, note that the mean radio-band SEDs of the 19 
sources for which L1216 ~ 10^"^ W Hz^^, which initially high- 
lighted a UV threshold, is no steeper than the SEDs of the 21-cm 
detections. This suggests that a bias towards larger radio sources is 
not the cause of the exclusive non-detections for the UV luminous 
sample. Also, presuming that the remaining five sources, for which 
the photon rate could not be determined, are below the threshold 
gives a detection rate (0 out of 5) which is not statistically signifi- 
cant, given the 40% chance of a detection at 2 < 1. 

Furthermore, at these redshifts we are probing look-back times 
more than double that of the 2: < 1 sources and if there is a 
larger fraction of smaller galaxies at these epochs, the gas will be 
ionised at lower luminosities. For instance, the highest photon rate 
for which there exists a 21-cm detection (5 x 10^^ s~^), is suf- 
ficient to fully ionise a gas distribution of scale-len gth of 1.6 kpc 
(50% t hat of a large spiral. 'Kalberl a & Kenjlioogl) . for no = 10 
cm ICurran & Whiting 2012)0 

If this is the case, objects even 
fainter than those targetted here are required to find the missing hy- 
drogen at high redshift, a task perhaps best suited for blind surveys 
with the Square Kilometre Array, where an optical redshift is not a 
prerequisite. 

To conclude, although we cannot rule out the radio structure 
and other hitherto unforeseen high redshift effects as the source of 
the non-detections, the fact remains that A ^ 912 A photons ionise 
hydrogen and there exists a 8.32 x 10~* probability (a 5.36cr sig- 
nificance) of the exclusive non- detections above a given ionising 
photon rate arising by chance jCurran & Whiting||2012l) . This is 
conjunction with the fact that the critical rate observed is just suf- 
ficient to ionise all of the gas in a large galaxy, leaves little doubt 
that photoionsation of the gas by the active nucleus is a dominant 
issue in the search of neutral gas within the hosts of these objects. 



For a disk, in which the thickn ess flares in the same manner as for 
the Milky Way kalberla et al]|2007l) . the total gas mass is estimated to be 
Mgas = 1.3 X 10^ Mq, a sixth of that foun d for a rate of 3 X 10^^ ion ising 
photons s^^, the critical value estimated bv lCurran & Whitind l2oT^ . 
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APPENDIX — SED FITTING 

In order to estimate the radio spectral indices and turnover fre- 
quencies, as well as the A ^ 912 A continuum luminosities 
and photon rates, we produced an SRD for each source by com- 
piling the NRD photometries, the Galaxy Rvolution Rxplorer 
(GALRX) near and far ultra-violet fluxes, as well as the BVRK 
magnitudes, published e lsewhere (compiled in lCurran et al. I I2OO8I : 
ICurran & Whitinsll2010t). After correc ting for Galactic extinction 
using the maps of lschlegel et al.l ( fl99^ , for each of the 241 sources 
the flux density was converted to a luminosity and the frequency 

de-redshifted to the source rest-frame^ 

Using Vandermonde matrices jPressetaLlI 19891) . three poly- 
nomials were fit to the data points for each source: 




10 15 

Rest— frame frequency, loQio^'^rest) ['~'^] 

Figure 7. SED fits to NVSSJ0205 10-1-224250 demonstrating a case where 
no radio turnover is apparent and where the A = 912 A continuum lumi- 
nosity cannot be determined (cf. Fig.|6). 

obtained from dLv /du = (with a maximum being tested using 
d^Ljy/df^ < 0) and the spectral index from the value of a = 
dL^/du at Urest = 1-42 x lO'' Hz, thus being defined as 5* oc u". 
Rig.[6]shows an example of a possible turnover (at 40 MHz) and a 
spectral index of a = —1.21. In the case where no radio turnover 
is apparent, providing dL,^/du is negative in the radio-band, the 
lowest observed frequency is assigned as an upper limit to this (as 
in Rig.lTj. That is, we assume that t he turnover occurs at a frequency 
lower than that currently observed jRanti et alJl99Clh . Otherwise, if 
dLv/dv ^ 0, the source is flagged as having an indeterminable 
turnover frequency. 

(ii) The full SED: A third order polynomial was fit to the all of 
the data points at j/rost ^ 10^' Hz. This fit provides a visual guide 
only. 

(iii) The UV SED: A first order polynomial (power-law) was 
fit to the 10"-^ < UrcBt < 10^^ Hz data. This upper cut-off is 
applie d since, although relevan t and used to derive the composite 
SRDs jCurran & Whitin3l2012h . the, usually sparsely sampled, X- 
ray points can have a strong influence on the fit, although these 
may be arising from other mechanisms than those producing the 
UV SRD: 

(a) Provided that dL^/dv < (e.g. as in Rig.|6j, the UV SRD 
was used to derive the luminosities and photon rate. 

(b) Should the gradient exceed dL^/dv = at Vrast = 
3.29 X 10^^ Hz, usually from a concentration of data points, 
the luminosity is estimated from the mean of any neighbouring 
data points (at v.^^t = iQis saio.so 

Only in the first case (a), where a reliable fit is obtained over the 
range specified, is the ionising photon rate derived from 



f 



■ dv, where logj^o Lv = ct logj^;, 1/ + C 



for a power-law fit, where a is the spectral index and C the inter- 
cept. This gives. 



lO'' 



lO'- 



-lO'^ 



where a < 0, 



(i) The radio SED: A third or second order polynomial was fit 
to the u-ccst ^ 10^^ Hz points. The turnover frequency was then 



h Jj^ ah " ah 

otherwise the ionising photon rate is flagged as indeterminable. 

Rinally, in Rig. [8] we show the SRD fits to the UV luminous 
sources (those which lie above the iuv = 10^'^ W Hz^^ cut-off. 
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Figure 8. As per Fig.|4] but for the 19 UV luminous sources 
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ICurran et alj|201ld) . 



the last row in Table [3]l. As seen from this, the ultra-violet lumi- 
nosi ties may be estimated ari d photon rates derived for all of these 
(see lCurran & Whiting|l2012h . 
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